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Abstract 



We determine the elements of the leptonic mixing matrix, without assuming 
unitarity, combining data from neutrino oscillation experiments and weak decays. 
To that end, we first develop a formalism for studying neutrino oscillations in 
vacuum and matter when the leptonic mixing matrix is not unitary. To be 
conservative, only three light neutrino species are considered, whose propagation 
is generically affected by non-unitary effects. Precision improvements within 
future facilities are discussed as well. 
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1 Introduction 



The experimental observation of neutrino oscillations both in the solar PQ-[1] and at- 
mospheric [5]- [8] regimes constitutes evidence for new physics beyond the Standard 
Model of particle physics (SM). Its favorite interpretation is neutrino masses, which 
allow flavour mixing. Neutrino masses and mixings have become an essential piece 
in the understanding of flavour dynamics. The determination of the absolute scale 
of neutrino masses, of its character -whether Dirac or Majorana-, as well as a deep 
and precise understanding of leptonic mixing and CP violation, are all complementary 
milestones required to unravel the flavour puzzle, which is one of the most important 
goals of elementary particle physics. 

Neutrinos are massless in the SM, since there are no right-handed neutrinos and 
the theory accidentally conserves (B — L), where B and L denote baryon and lepton 
number, respectively. Neutrino masses can be Dirac or Majorana. Pure Dirac neutrino 
masses added to the SM are disfavored theoretically for the following reasons: 1) They 
worsen the problem of the large mass hierarchy between neutrinos and other fermions, 
with Yukawa couplings down to 11 orders of magnitude smaller for the former. 2) They 
are unnatural in the 't Hooft sense, as Majorana mass terms of right-handed neutrinos 
are not forbidden by the SM gauge symmetry; Dirac neutrinos require to impose the ad 
hoc assumption that the Lagrangian conserves L, and consequently (B — L) so long as B 
is conserved. 3) Such conservation precludes the simplest scenarios for the generation 
of the baryon asymmetry of the universe through leptogenesis [H], at energies above 
the electroweak scale. 

Non-zero neutrino masses are thus a signal of physics beyond the SM. The new 
scale associated to it, M, is likely to be much higher than the electroweak scale. When 
M ^> Mz, the effects of new physics at the high energy scale can be parameterized at 
low energies, without loss of generality, by an effective Lagrangian including: 

• Corrections to the parameters of the SM Lagrangian. 

• The addition to the SM Lagrangian of a tower of non-renormalizable -higher 
dimensional- effective operators, which are SU(3) x 577(2) x £7(1) gauge invariant. 
Their coefficients are suppressed by inverse powers of the large energy scale M. 

The relative intensity of the different higher- dimensional effective operators is model 
dependent. The SM couplings themselves are expected to be generically modified, 
though, in any scenario involving new physics at high energy. We concentrate on 
modifications to those couplings, in the present work. 

Among the SM couplings, the PMNS mixing matrix is introduced in connection 
with charged current interactions of leptons. Its departure from the unit matrix is the 
origin of all leptonic mixing and putative CP- violation effects in neutrino physics, which 
are essential building blocks of the flavour puzzle. In typical analyses of experimental 
data, the mixing matrix is assumed to be unitary. 
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The complete Lagrangian including new physics is necessarily unitary and proba- 
bility is conserved. Unitarity violation is one typical low-energy signal of models of 
new physics, though, when data are analyzed in the framework of the SM gauge group 
with three light species of quarks, leptons and neutrinos, plus mass terms for the latter. 
Notice that for the analogous matrix in the quark sector, the CKM matrix, deviations 
from unitarity are considered a good window for physics beyond the SM and extensively 
studied. 

An evident possible source of non-unitary effects in leptonic mixing is the hypothet- 
ical existence of more than three light neutrino species, that is, light sterile [10] neutri- 
nos. Their theoretical implementation typically requires strong fine-tunings, though, 
and no definite signal exists for their existence [H]-[I3]. To be conservative, we will 
disregard such possibility in the analysis below. Would light sterile neutrinos turn out 
to be indeed present in nature and detected, its presence in mixing processes would 
obviously lead to stronger signals than those developed in this paper. 

Even in the framework of only three light neutrinos, extensions of the SM postulated 
in order to generate the observed small neutrino masses typically produce a leptonic 
mixing matrix which is non-unitary. For instance, any theory in which the neutrino 
mass matrix turns out to be part of a larger matrix, which involves heavy fields, may 
generically result in an effective low-energy non-unitary mixing matrix for the light 
leptonic fields [E]-[I7]- A paradigmatic example is the seesaw mechanism [18]. In it, 
light neutrino masses take natural (i.e. non-fine-tuned) values, in contrast to all other 
fermion masses except the top quark mass. Heavy sterile (right-handed) neutrinos, 
with masses far above the electroweak scale, are introduced at the new high energy 
scale. Although the non-unitary mixing induced at low energies in the canonical -type 
I- seesaw model is typically expected to be too small for detection, this is not necessarily 
true for variants of the seesaw mechanism, or other theories beyond the SM. 

At the same time, we are about to enter an era of high precision neutrino physics. 
With on-going and forthcoming experiments, as well as the future facilities under dis- 
cussion, future neutrino oscillation experiments will be aiming at a measurement of 
the last unknown leptonic mixing angle #13, even if very small, as well as of leptonic 
CP violation. It is pertinent to ask whether such precision can shed further light on 
unitarity. Departures from it would probe the new physics behind. 

We will thus relax here the assumption of unitarity in the low-energy leptonic 
mixing matrix of weak interactions, and let data rather freely tell us up to what point 
the measured elements of the mixing matrix arrange themselves in a unitary pattern. 
This will allow to identify the less constrained windows in flavour space and thus the 
most sensitive ones, as regards new physics. 

We will not work in any concrete model of neutrino masses. Nevertheless, as only 
neutrino masses clearly signal new physics -in contrast to masses for charged leptons 
or the rest of the fermions-, it is plausible that the new physics behind sneaks through 
at low energies primarily through its effects on neutrino propagation. That is, in this 
first work we implicitly assume that the physics of fields other than neutrinos will be 
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that of the SM. We can summarize then our approximations on a set-up that we will 
dub Minimal Unitarity Violation (MUV), based on the following assumptions: 



• Sources of non-unitarity are allowed in those terms in the SM Lagrangian which 
involve neutrinos. 

• Only three light neutrino species are considered. 

Leptonic and semileptonic decays, together with neutrino oscillations, will be analyzed 
in this minimal set-up. Supplementary non-unitary contributions to physical tran- 
sitions can result from new physics affecting the SM couplings, for fields other than 
neutrinos and/or higher- dimensional operators in the effective Lagrangian. Barring ex- 
treme fine-tunings, they should not affect the order of magnitude of the results obtained 
in the MUV scheme, though. 

The paper is organized as follows: Sec. [2] defines the framework and introduces the 
non-unitary mixing matrix N which replaces the unitary PMNS matrix. In Sec. [3l a 
formalism is developed for the study of neutrino oscillations in vacuum and matter, with 
non-unitary leptonic mixing. Sees. H] and [5] deal with data re-analyzed in the MUV 
scheme: present data from neutrino oscillation experiments are considered in Sec. H] 
and the mixing matrix resulting from their analysis is obtained, while the unitarity 
constraints resulting from W and Z-decay data, lepton universality tests and rare 
charged lepton decays are presented in Sec. [51 The final mixing matrix resulting from 
the combination of oscillation and weak decays data is presented in Sec. |6j In Sec. 
the impact of future experiments on the results obtained in the previous sections is 
studied. In Sec. El we conclude. 



2 Non-unitary mixing matrix 

Our aim is to substitute the usual leptonic unitary matrix Upmns by a non-unitary one. 
To be definite, we will analyze the consequences of an effective low-energy Lagrangian 
which, in the mass basis, reads 



C eff = ~ {vii^Vi - v c (mi Vi + h.c.) - ^= (W+ l a 7 M (1 - 75) N ai v { + h.c.) 
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(Z„ 9i Y (1 - 75) (rtN)ij + h.c.) + ... 



2 cos 9w 

where z/j denotes four- component left-handed fields. Eq. ([I]) is the usual Lagrangian 
for neutrinos in the mass basis, albeit with the unitary matrix Upmns m the charged 
current substituted by a general non-unitary matrix N. The neutral current coupling 
has been modified as well, as expected in general from the non-unitarity of N. Notice 
that the Lagrangian includes a Majorana mass term for neutrinos, for the sake of 



3 



definiteness, although for our numerical analysis below it would make no difference to 
consider neutrinos of the Dirac type. 

It is necessary to clarify the relation between mass and flavour eigenstates. With 
non-unitarity present, both the mass and flavour bases cannot be orthonormal. Indeed, 
N connects the quantum fields in the mass basis with those in the flavour basis3 where 
the weak couplings are diagonal, 

v a = N ai vi . (2) 
While the canonical kinetic terms in Eq. ([T]) give rise to orthonormal mass eigenstates, 

(vi\vj) = 8ij> (3) 
consistency between quantum states and fields requires the identification [T9] 

^--7wm^ KM) ^ :M) ' (4) 

where, on the right-hand side, the normalization factor has been absorbed in the def- 
inition of N. This equation is to be compared with Eq. (j2j). It follows from it that 
flavour eigenstates are no more orthogonal^: 

{vpW) = {N^)p a ^8 aP , (5) 

which will induce relevant physical effects, as it will be shown later on. 

Prior to any predictions for physical transitions, the Lagrangian must be renor- 
malized, determining from experiment its free parameters. In particular, the weak 
coupling in Eq. ([T]) differs from the SM expression. Accordingly, the Fermi constant 
measured in experiments cannot be identified anymore with the SM tree level combi- 
nation Gf = V / 25 ,2 /(8M^ / ), due to non-unitarity. For instance, the Fermi constant Gp 
extracted from the decay fi — > u^ePe is related to Gp by 

G F = 7=^= ==- (6) 



y/(NW) ee (NW) 



ft' 



The rest of the parameters of the Lagrangian coincide with those in the standard 
treatment. 



throughout the paper, Greek (Latin) indices label the flavour (mass) basis. 

2 Notice that these are effective low-energy flavour eigenstates. In the corresponding complete -high 
energy- hypothetical theory, it should be possible to define an orthonormal flavour basis. 
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The effective Lagrangian in the flavour basis 

Let us write the Lagrangian in Eq. ([1]) in the flavour basis. Upon the transformation 
in Eq. (j2J), it follows that 

£ eff = ~ (iv a HNN^)~l v p - ^ Q [{N^fmN-XpVp + h.c.) 
-^=KU"(l-7> a U.c.) 

- 5-^- (^P a y (1-75)1/0+ h.c) + ... , (7) 

Z COS u\y 

where m = diag(mi, 7712, 7713). In this basis, weak couplings are diagonal and a general 
neutrino mass matrix and kinetic term have appeared. A neutrino mass term can 
always be diagonalised by a unitary transformation (or bi-unitary in the case of Dirac 
neutrinos), leading to a unitary contribution to the mixing matrix, as in the standard 
treatment. Kinetic terms must be both diagonalized and normalized, though, so as to 
obtain canonical kinetic energies. 

The hypothetical different normalizations induced in the kinetic energy of neutrino 
fields by the new physics, is the key point to obtain non-unitary effects in neutrino 
mixing, in the MUV scheme. Whenever at least two normalizations of neutrino fields 
differ, a non-unitary weak mixing matrix follows. 

For the sake of physics intuition, let us discuss the theoretical consistency and 
implications behind the Lagrangian in Eqs. ([I]) and (jTJ), although the hypothetical 
reader mainly interested in the numerical analysis is invited to proceed directly to the 
next Section. 

A pertinent question is whether there exists a 577(2) x £7(1) invariant formulation 
of the MUV formulation above, as it should. Consider a generic effective Lagrangian 
valid at energies less than a high scale M of new physics, M 3> Mz, resulting after 
integrating out the heavy fields present above such scale. The effective Lagrangian has 
a power series expansion in 1/M of the form 

C*ff = C SM + SC d=5 + 5C d=6 + ... , (8) 

where Csm contains all S77(3) x SU{2) x U{1) invariant operators of dimension d < 4 
and the gauge invariant operators of d > 4, constructed from the SM fields, account 
for the physics effects of the heavy fields at energies < M. After electroweak sponta- 
neous symmetry breaking (SSB), the operators with d > 4 will give corrections to the 
couplings present in the SM Lagrangian and also produce new exotic couplings. 

The question is whether there exist 577(2) x f7(l) gauge invariant formulations of 
the Lagrangian which give rise, after SSB, to flavour mixing corrections to the neutrino 
kinetic energy, as discussed here. This is the case indeed. There exists for instance 
a d — 6 gauge invariant operator which precisely results in the corrections in the 
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MUV scheme. It is the operator characteristic of the canonical seesaw mo del and its 
generalizations [20], and also of some extra-dimensional constructions [21] . 

6C d=e = c d a f (Lj)ip(jpL p ), (9) 

where L denotes left-handed leptonic doublets^, c d=6 is the -model dependent- coeffi- 
cient matrix of 0(1/ M 2 ) and is related to the standard Higgs doublet <fi by <fi = ZT20*. 

Other SU(2) x U(l) invariant operators may be written as well, resulting generically 
-after SSB- in corrections to both the neutrino and charged leptons kinetic energy. 
There exist even SU(2) x U(l) invariant operators contributing to the charged lepton 
kinetic energy and not to that of neutrinos, such as 

6C d = 6 = c' a f 6 (L a <f>) (tfLp) , (10) 

with c /d=6 being the coefficient matrix of order 0(1/M 2 ). Theories with Yukawa cou- 
plings to heavy extra fermions -be it of Dirac or Majorana type- can easily give rise 
to such effective couplings at low energy. Such an operator leads -after SSB- to a La- 
grangian with the same couplings to the W boson as in Eq. ([I]) , albeit with corrections 
in the Z-charged lepton couplings instead of in the Z-neutrino ones. It means that 
all results obtained below from W exchange alone would also hold for the purpose of 
constraining such theories. 

Many other gauge invariant operators [22] can contribute to the MUV physics 
discussed here. The last comments above are only a digression, to illustrate the rather 
general application realm of our treatment, although we will keep our focus on neutrinos 
as stated. 

In summary, a generic model is expected to give rise after SSB to modifications of 
the standard couplings, as well as new exotic ones. In the minimal scheme analyzed 
in this paper, MUV, only the former will be taken into account and more precisely 
only those couplings involving neutrinos as specified above, as it is natural to expect 
that the latter will be specially sensitive to the new physics responsible for neutrino 
masses. This simplification should provide a sensible estimation of the best windows for 
non-unitarity, unless extreme fine-tunings and cancellations occur between the different 
type of contributions, in some hypothetical model. 



3 Neutrino oscillations without unitarity 

Let us consider now the impact of the Lagrangian in Eq. ([T]) on neutrino oscillations, 
both in vacuum and in matter. 

3 In addition to the well-known d = 5 operator responsible for neutrino masses, SC d=5 = 
I c % 5 {paft^j (ft Lp\ + h.c. , where c^ 5 is the coefficient matrix of 0(1/M). 

4 As the factors in parenthesis in this equation are singlets of SU(2) x U(l), ft is tantamount to p 
in this operator. 
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3.1 Vacuum oscillations 



Consider free neutrino propagation, described by the free Hamiltonian H^ ree , resulting 
from the first two terms in the Lagrangian, Eq. (pQ). The time evolution of mass 
eigenstates follows the usual pattern. Indeed, 

= (11) 

and because of the orthogonality of the mass basis, 

(vjlH^M^SijEi, (12) 

where are the eigenvalues. Using now the completeness relation, J2j \ u j)( v j\ = 1) 
Eq. dHJ reads: 

l J t \^ = E \vj)(vj\H free \vi) = EM) , (13) 

3 

which is the usual time propagation for free states. 

Consider now instead the free evolution in the flavour basis, which is not orthonor- 
mal and for which there is not the usual completeness relation, as ^ a \v a )(v a \ ^ 1. 
The time evolution is given by 

4k) = # /ree k>, (i4) 

at 

which, using the orthogonality and completeness of the mass basis, results into 
d 



% Jt' y ' 



'«> = £ WM^H^Wa) = ^(iV* E (NT'UW?) > (15) 
3 P 



where E = diag(£?i, E2, E%). This is to be compared with the A^-dependence of the 
matrix elements between flavour eigenstates, given by 

(vp\Hf™\v a ) = {N*EN l U. (16) 

That is, the evolution in flavour space is driven by the combination (N* E (iV*) _1 ) 
and not by the product (iV* EN 1 ) appearing in Eq. (|16p . in contrast to the customary 
expression in standard -unitary- treatments. Because of the non-unitarity of iV both 
expressions are no more equivalent. Technically, this is a key point in the different 
results for the non-standard case, to be obtained below. 

Notice, moreover, that the combination (iV* E (iV*) -1 ) is not Hermitian, even if 
the free Hamiltonian itself is Hermitian. This in turn implies that the evolution of 
flavour bra states, (u a \, differs from the evolution of the flavour kets, leading both to 
the same probability equation, as they should. 
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The analysis of free propagation in space is analogous to that for time evolution 
described above and we will not repeat it in detail. Flavour eigenstates, after a distance 
L, transform into 

\v a {L)) = ^N* ai e i » L (N*)J\vJ, (17) 

if 



where Pj are the momentum eigenvalues, Pi — yEf — mf. The oscillation probability 
after traveling a distance L can now be obtained, 

I IV* piPiL M 1 2 

P VaV0[ E,L) s = %^ aiNN Z = W 

+2Re {N al N^N^Np2} cos A 12 + 21m {N al N^N^N^} sin A 12 + 
+ 2Re {N a2 N* 3 N; 2 Nf) 3 } cos A 23 + 21m {N a2 N* 3 N; 2 Np 3 } sin A 23 + 
+ 2Re {N^N^N^Npx} cos A 31 + 21m {N a3 N* al N; 3 N pi } sin A 31 

where Ay = AmfjL/2E, with Am| = — m|, as usual. Written in this way, the 
expression is easily seen to reduce to the standard one if N was unitary, as it should. 

The first very important consequence of Eq. ffl8l) is that the non-unitarity of iV is 
shown to generate a "zero-distance" effect [15], i.e. a flavour transition already at the 
source before oscillations can take place. Indeed, for L = 0, it follows that 



|(AWt) 
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p - |s ^^»u^' (19) 

an effect that can be tested in near detectors, thus setting strong limits on unitarity 
as we will see later. Nevertheless, due to non-unitarity, the probability as defined in 
Eq. (JT81 does not sum up to a total probability of 100%. To make contact with data, 
let us discuss the implications of our treatment for the production and detection cross 
sections and, finally, for the number of events detected in a given experiment. 



3.1.1 Production/detection cross sections and widths 

The non-unitarity of the mixing matrix N implies the following corrections, for pro- 
cesses computed at tree-level: 

• Charged current (CC) cross-sections and fluxes involving only one neutrino flavour 
a are given by 

a cc = a cc { SM) {NN%a = {NN%a ( (20 ) 
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where <Ja C ^ SM ^ and ^ c<ySMS) are the SM cross section and flux, respectively. The 
same correction factor affects decay widths involving one neutrino flavour. 

Charged current cross sections involving two neutrino flavours, a, (3, will be mod- 
ified into 

*2? = <? 5M) (NN% a (NN% P , (21) 

with the same weight factor affecting widths or fluxes involving two neutrino 
flavours. 

Neutral current (NC) processes are weighted by a different combination. A decay 
width involving two neutrino mass eigenstates, Uj, is given by 

r(Z -> v iVj ) = T SM (Z - v iVi ) I (JVtiV)y | 2 . (22) 

Analogously, when detecting a neutrino z/j through neutral current interactions, 
as in SNO, modified cross sections will have to be considered, 



^<7^ 5M )|(ArtA%| 2 , (23) 



where the sum over j is due to the fact that the final neutrino Uj remains unde- 
tected. 



3.1.2 Number of events 

The number of events in a detector located at a distance L away from the source 
would be given, apart from backgrounds, by the convolution of the production flux, the 
oscillation probability, the detection cross-section and the detector efficiency, integrated 
over energy. In short, 

n event s ~ J dE Pu a u (E, L) a p {E) e(E), (24) 

where d$ a (E)/dE is the neutrino flux, <Jp(E) is the detection cross section and e(E) 
the detection efficiency In the presence of MUV, all factors in Eq. (^/^) should be 
corrected, as discussed above. It is easy to see that there are cancellations between the 
different iV dependent factors they exhibit. 

For instance, for experiments in which both production and detection take place 
via charged currents, involving each one neutrino flavour, the denominator of P VaVf) 
-Eq. (1181) - cancels the correction factors in the flux and cross section, Eq. ( |20|) . This 
allows to express in this case the number of events simply as 

n events ~ J dE dE K ' Ku, (L, E) a C p C{SM) (E) e(E) , (25) 
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where P UaU0 (L, E) is the probability in Eq. f|T8|) . amputated from the normalization 
factors in its denominator, 

P UaU0 (L,E) = \J2N* m e lP * L N^\ 2 . (26) 

i 

It turns out that, in practice, most experiments extract the probabilities from the 
measured number of events, parameterized -via Monte Carlo simulations- in terms of 
the SM fluxes and cross sections, precisely as in Eq. (I25I) . Within MUV, their analysis 
thus provides a direct estimation of P UaU/3 (L, E) in Eq. (126]) . This is the case for the 
very large number of experiments in which neutrinos are detected via charged current 
interactions and produced from decays of hadrons like 7T, K or (3 decays. 

Obviously, there are exceptions. For instance, if the neutrino flux expected in the 
far detector of the previous example is not taken from a Monte Carlo simulation, but 
from a direct measurement in a near detector, the cancellation described in the previous 
paragraph would not be complete and extra iV-dependent factors will have to be taken 
into account. 

Besides, when the production mechanism is not hadronic, but leptonic, as from 
/i or even r decays, the fluxes would need two corrections instead of one, since the 
production involves two insertions of N, as in Eq. (12"T1) . For instance, the neutrino 
fluxes produced at a Neutrino Factory from /i decay should thus be corrected by the 
factor (NN^^(NN^) ee . 

Finally, the analysis of detection through neutral current processes is modified as 
well. Since such processes are sensitive to the sum of all neutrino species, the number 
of events is given by 

/ ifaOO\SM) / t-i\ 
i 

with (E) as in Eq. (|23l) and -for instance for propagation in vacuum- 

P VaVi {L,E) = \N m \ 2 . (28) 



3.2 Matter effects 

When neutrinos pass through matter, they interact with the medium and the effect is a 
modification of the evolution and consequently of the oscillation probability. Consider 
the usual interaction Lagrangian 

- C mt = v / 2G F n e ^ e 7 z/ e ^=G F n n ^ z/ Q 7o^ , (29) 

where n e and n n are respectively the electron and neutron densities. The first term 
corresponds to charged interactions, while the second term corresponds to neutral 
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interactions. In the mass basis within MUV, Eq. fl29l) reads 

- C mt = V cc N *eiN ej moVj - V NC KiN^H^ , (30) 



where Vcc — V2Gpn e and V^c = -^Gpn n . In order to estimate the time evo- 
lution of states passing through matter, consider the interaction Hamiltonian, H mt , 
corresponding to Eqs. (I29p and (I30p . Its matrix elements in the mass basis read 

Ef = {vj\H int \vi) = V cc N ei N* ei - Vnc{N^ N)ji , (31) 



or, m matrix notation 



H mt = [iVt diag(V cc - V NC , -V NC , -V NC ) N] 1 . (32) 
The evolution equation for mass eigenstates in matter is then given by 

. d 



where E is the energy matrix for free eigenstates, introduced in Eq. ( TT3|) . In contrast, 
the evolution through matter of flavour eigenstates is given by 



d , 
l dt lUc 



) = J2[N*(E + H mt )(N*)- 1 ] Wis), (34) 



where again (N*) 1 cannot be traded by N f -as it is usually done in the standard case- 
because N is not unitary. 

It is easy, although cumbersome, to write now explicitly the equations above for 
the three family case. To illustrate the main new effects, it is enough to write here 
explicitly the effective flavour potential in the second term in Eq. (j34j) . for the case of 
two families: 

ft*H int (N*)- 1 = N*N t ( Vcc ~ VNC ° )jV*(iY*)- 1= (35) 

V U —VNC J 



(Vcc - V NC )(NN% e -V NC yJ$g}g{NNl) 



fie 



Consequently, MUV results generically in exotic couplings in the evolution through 
matter. This effective potential is not diagonal, in contrast to the unitary case. More- 
over, the neutral current contribution can not be rewritten as a global phase in the 
evolution equation and thus it contributes to the oscillation probabilities. 
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4 Matrix elements from neutrino oscillations 



We will use now the most relevant data on neutrino oscillations, to determine the 
elements of the mixing matrix, without assuming unitarity. In this first work, we do 
not perform an exhaustive analysis of all existing oscillation data; our aim is rather 
to estimate what is the role played by the different experiments in constraining the 
matrix elements. 

All positive oscillation signals available nowadays^] correspond to disappearance 
experiments [Tj-[S]. Since the disappearance oscillation probability in vacuum is given 
by 

P VaVa = |iV Ql | 4 +|iV a2 | 4 + |iV Q 3| 4 + 2|iV al | 2 |iV a2 | 2 cosA 12 

+ 2|iV Ql | 2 |iV Q 3| 2 cosA 1 3 + 2|iV Q2 | 2 |iV a 3| 2 cosA 2 3, (36) 

disappearance experiments may provide information on the moduli of elements, while 
phases will remain unknown, as in the unitary case. Furthermore, as no v T disappear- 
ance experiment has been performed yet, this type of vacuum experiments will only 
constrain the elements of the e and the //-rows, as Eq. (|36|) indicates. Nevertheless, 
the no-oscillation results from some appearance experiments will also provide useful 
non-unitarity constraints. 



Vacuum oscillations 



The exact appearance and disappearance probabilities in vacuum, Eqs. (1TB]) and (1261) . 
will be used in the numerical analysis. To illustrate the discussion, the amputated 
probabilities can be approximated as follows, though, for some experiments studied 
below, depending on the range of L/E: 



A12 — 0. Eq. fl36l) reduces then to -for instance for the case of v e disappearance 
and disappearance- 

P PePe ~ (|Ar el | 2 + |iV e2 | 2 ) 2 + |iV e 3| 4 

+ 2(|iV el | 2 + |iV e2 | 2 )|iV e 3| 2 cosA 2 3, (37) 



~ (|iv Ml | 2 + |iv M2 | 2 ) 2 + |iv M 3| 4 

+ 2(|iV Ml | 2 + |iV M2 | 2 )|iV / ,3| 2 cosA 2 3, (38) 

respectively. Relevant experiments in this class include CHOOZ [1], a reactor 
experiment sensitive to v e disappearance, as well as the disappearance atmo- 
spheric [6] and accelerator experiments such as K2K [7] or MINOS [8j. Eqs. fl37|) 

5 Except for the LSND experiment [11] which is currently being tested by MiniBooNe [13] and 
which we will not consider in this paper. 
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and (|38|) indicate that this type of vacuum experiments cannot disentangle by 
themselves the element \N a i\ from l-A^I, as they appear in the combinations 
|iV el | 2 + |iV e 2| 2 and | ^V Ml | 2 + \N^ 2 \ 2 , respectively. In addition, the equations show 
as well the presence of a degeneracy between those combinations versus |A^ e3 | and 
|A^ 3 |, respectively 

• A12 7^ with A23 3> 1. The latter -atmospheric- oscillation frequency is aver- 
aged out resulting in -for instance for v e disappearance- 

P PePe ~ |iV el | 4 + |iV e 2| 4 + |iV e 3| 4 + 2|iV el | 2 |iV e2 | 2 cosA 12 . (39) 

KamLAND [3] is a reactor experiment with a longer baseline than CHOOZ and 
falling into this category. Notice that the dependence on |iV el | and |-/V e2 | in 
Eq. f)39p differs from that in Eq. (13 7p . suggesting that the combination of both 
type of experiments may help to tell those elements apart, as it will be shown 
later on. 

• A 12 — and A 2 3 — 0. The appearance and disappearance probabilities corre- 
spond then to a simple formula (see Eq. (fT9l ): 

P VaV0 ~ KiVivt)^! 2 . (40) 

KARMEN [12] and NOMAD [23] are appearance experiments in this class, well 
described by Eq. (TJUJ); the same holds for the data on disappearance at the 
near detector in MINOS and on v e disappearance at BUGEY [21]. 

Oscillations in matter 

A very important experiment in this class is SNO. In the unitary treatment, the v e 
produced at the core of the sun are approximately eigenstates of the total Hamilto- 
nian, since the interaction Hamiltonian dominates the evolution in this region of the 
sun. Denoting the total eigenstates in matter by |z>j), a v e in the center of the sun is 
approximately [2"5] : 

\v e ) ~ v / 0T|£i> + y/0H\v 2 ), (41) 

within a 2% accuracy. The state then evolves adiabatically so that, when leaving the 
sun, the \vi) states can be replaced by the vacuum eigenstates |z/j), leading to 

P VeVe ~ 0.1|iV el | 2 + 0.9|iV e2 | 2 , (42) 

which allows a clean measurement of |iV e2 | 2 . 

Within the MUV scheme, a priori the analysis varies. This was illustrated in 
Eq. (135]) for two-family oscillations in matter, which exhibits exotic non-diagonal terms 
and where the neutral currents may play a priori a significant role. Nevertheless, we will 
see below that the absence of oscillation signals at near detectors constrain deviations 
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from unitarity, for all (NN^) elements but (NN^) TT , to be smaller than O(10~ l ). In 
fact, it will turn out that all bounds on (NN^), including (NN^) TT , are improved 
also from weak decays and the values of the off-diagonal elements constrained to be 
smaller than 5%, as it will be shown in the next Section. In consequence, for the level 
of precision aimed at in this work, it is unnecessary to perform the complete MUV 
analysis of SNO data and Eq. (142]) keeps being an appropriate approximation. This 
determination of | iV e 2 1 will be a major input in resolving the MUV degeneracy between 
|JVei| and \N e2 \. 

In all numerical analysis below, the values of Ami 2 and Am 2 3, resulting from our 
fits in the MUV scheme, will not be shown: they coincide with those obtained in the 
unitary treatment, as expected from the fact that the oscillation frequencies are not 
modified in the MUV scheme, unlike the amplitudes. 

4.1 Constraints on the e-row 

In Fig. [T] (left) we present the 1, 2 and 3a contours of a three-family fit to CHOOZ 
data, combined with the information on Am^ resulting from an analysis of K2K data. 
The dotted line represents the unitarity condition (NN^) ee = 1. 

Since CHOOZ data are compatible with the no-oscillation hypothesis, the fit shows 
allowed regions in which the first line in Eq. (1371) is close to one, while the second 
-oscillatory- term vanishes. That is, either |iV el | 2 + \N e2 \ 2 ~ 1 with |iV e 3| 2 ~ 0, or 
\N el \ 2 +\N e2 \ 2 ~ with \N e3 \ 2 ~ 1. The detection of the L/E dependence in KamLAND 
selects the first combination, though, see Eq. (1391 . The significant loss of sensitivity to 
|iV e 3 1 of the 3a contour with respect to the 1 and 2a ones can be understood from the 
fact that CHOOZ loses its sensitivity for Am| 3 ~ 0.001, as can be seen in Fig. 55 of 
Ref. [1]. Indeed K2K excludes such small values of Am^ at 1 and 2a, but not at 3a, 
where the loss of sensitivity occurs. Notice that the 3a contour intersects the unitarity 
condition at |iV" e3 | 2 ~ 0.05, which agrees with the usual bounds for |iV" e3 | 2 , obtained 
under the assumption of unitarity. 

KamLAND also helps to disentangle |iV e i| 2 from 1 7V e 2 1 2 , as illustrated in Fig. [1] 
(right), in which the 3a contour of a fit to KamLAND data is presented (dashed line), 
combined with those from CHOOZ and K2K. Since CHOOZ constrained | iV e i | 2 H- 1 JV e 2 1 2 
to be close to 1, only a narrow strip near the diagonal is allowed. The region is still 
large, though, due to the symmetry of Eq. (1311 under the interchange of |iV el | 2 with 
|iV e2 | 2 . 

This final degeneracy can be lifted with information from SNO. The SNO data 
on the ratio of the charged-current over neutral-current fluxes results in the rather 
horizontal 3a strip (dotted line) in Fig. [1] (right). To determine this region, the ratio 
of charged-current over neutral-current fluxes [2] can be approximated by Eq. (j4"2l . A 
5% variation has been allowed, to take into account the corrections stemming from 
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Figure 1: Left: 1, 2 and 3a contours for CHOOZ and K2K data (solid lines) with 
unitarity condition (dotted line). Right: 3a contours for CHOOZ, K2K and KamLAND 
data (dashed line), SNO data (dotted line) and their combination at 1, 2 and 3a (solid 
lines). 



Eqs. ( 1271) and (|35l) . Furthermore, we have verified that even a 10% variation in those 
coefficients would not change significantly the results of the fit. 

The combined fit of CHOOZ, KamLAND, SNO and K2K data is depicted at 1, 
2 and 3a by the solid contours in Fig. [1] (right). The Figure shows then that the 
combination of all this complementary information constrains all elements of the e-row 
with a precision only slightly inferior to that of the usual unitary analysis. 



4.2 Constraints on the /Li-row 

In Fig. [2] we show the 1, 2 and 3a contours (solid lines) of a fit to K2K data, combined 
with an estimation for SK. The latter resulted from translating the measured value of 
sin 2 (2^23) in Ref. [6] to matrix elements, using Eq. fl38|) . The dotted line represents the 
unitarity condition, (NN^)^ = 1. Without additional information at different L/E, 
liV^xl 2 and |A^ 2 | 2 can not be disentangled and our knowledge of the /i-row is much 
worse than when imposing unitarity 

Putting together all the information developed above from different oscillation ex- 
periments, the following allowed ranges are obtained -at 3a-, for the elements of the 
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Figure 2: 1, 2 and 3a contours for K2K and SK data (solid line) with unitarity condi- 
tion (dotted line). 



leptonic mixing matrix: 

0.75 - 0.89 0.45 - 0.66 < 0.34 

\N\= \ [(|iV M i| 2 + |iV Aia | 2 ) 1 /2= 0.57- 0.86] 0.57-0.86 ]. (43) 

? ? ? 

We observe that, without the assumption of unitarity, present oscillation signals 
can only determine half of the matrix elements. The elements of the first row have 
all been determined or constrained (|A e 2| 2 mainly by SNO, (A^l 2 mainly by CHOOZ 
and |iV e i| 2 by KamLAND combined with the others). In contrast, for the second row, 
present atmospheric and accelerator experiments are unable to discriminate between 
jA^x) 2 and |A^ 2 | 2 - Finally, no direct information can be obtained on the r-row in the 
absence of u T oscillations signals. 



4.3 Constraints on non- unitarity from near detectors 

NOMAD, KARMEN, BUGEY and the near detector at MINOS provide constraints 
on the elements of NN^ as it follows from Eq. (|4*0|) . We obtain, at the 90%CL, 

'1.00 ±0.04 <0.05 
|AW f |Ri[ < 0.05 1.00 ±0.05 < 0.013 | . (44) 
< 0.09 < 0.013 




With this information, |A^i| and (A^l in Eq. (I43p can now be disentangled. All 
in all, the mixing matrix resulting from analyzing oscillation data within the MUV 
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scheme is given by 

/ 0.75 - 0.89 0.45 - 0.66 < 0.27 \ 

\N\=\ 0.00 -0.69 0.22- 0.81 0.57- 0.85 . (45) 

\ ? ? ? / 

Notice that, even adding the constraints obtained at near detectors, not all matrix 
elements can be determined from oscillation data. 



5 Constraints on non-unitarity from electroweak 
decays 

Neutrino oscillations are evidence of a non-trivial leptonic mixing, allowing to deter- 
mine the individual elements of the mixing matrix from its data, as done in the previous 
Section. In contrast, leptonic and semileptonic decay data are not appropriate for this 
task. This is because, contrary to the quark sector, where the different quark mass 
eigenstates can be tagged, neutrino eigenstates are not detected separately. The exper- 
imentally measured rates correspond then to sums over all possible mass eigenstates, 
resulting only in sums of products of matrix elements. 

Leptonic and semileptonic decays may be sensitive to leptonic non-unitarity, though, 
due to the "zero-distance" effect, which in the flavour basis is encoded by its non- 
orthogonality, Eq. (jSJ). The combinations (NN<) a p can be extracted from them, as 
suggested by Eqs. f[20l) - fl23l) . With that aim, W, Z, ir and lepton decays are analyzed 
in this chapter, in the MUV scheme. The results will further constrain the mixing 
matrix obtained from neutrino oscillation processes. 



5.1 W decays 

With a non-unitary leptonic mixing matrix N, the decay widths for W into charged 
leptons and neutrinos are given -as in Eq. (|20|) - by 

r( W - t a v a ) =Vr(f^ e a vi) = -^=^(NN% a . (46) 

. 6v27r 

Gp has been related to the Fermi constant Gp, measured from the decay \\i — > v^ev^ 
by Eq. ([6]), allowing to extract now from Eq. (|46p the following combinations: 

(NN% a T(W -> i a u a ) 6V2tt ^ f 
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M 



Using the results for the W decay widths and mass from Ref. [26J, as well as G F 
(1.16637 ± 0.00001) ■ 10" 5 , the parameters f a are 

f e = 1.000 ±0.024, 
/„ = 0.986 ±0.028, 

f T = 1.002 ±0.032. (4* 



5.2 Invisible Z decay 

Further constraints stem from the invisible Z-decay width, which, for non-unitary 
leptonic mixing N, is given by (see Eq. (1231 ) 

G A/f 3 

V{Z - invisible) = £ V{Z -> v iVj ) = ^jJ- (1 + Pt) £ > ( 49 ) 

where pt ~ 0.008 [26J takes into account radiative corrections mainly steming from loops 
with the top quark. As the dominant radiative corrections do not involve neutrinos, the 
dependence on the mixing matrix in Eq. (j49p appears as a global factor to an excellent 
approximation. Using Eq. (jHJ), the equality ^ • \(N^N)ij\ 2 = Y^a,p \{NN^) a p\ 2 and 
the data provided in Ref. [26], the following constraint is obtained 

E^Kaw'M 2 _i2V2.r ( z^ invisible) = 2 . 984±a009 . (50) 



y/{NW)„(NM) m GjfM&l + ft) 

As it is well known, this number should correspond to the number of active neutrinos 
at LEP. Its 2a departure from the value of 3 is not (yet) significant enough to be 
interpreted ClS db SI gnal of new physics @. 



5.3 Universality tests 

In addition, ratios of lepton, W and 7r decays, used often in order to test universality [261 
[27] , can be interpreted as tests of lepton mixing unitarity. They result in constraints 
for the diagonal elements of NN\ as suggested by Eqs. f[2"0l - f[2"3l and resumed in Table 
I. 

The processes investigated so far constrained the diagonal elements of the product 
NN^. Limit values for its off-diagonal elements can be obtained instead from rare 
decays of charged leptons, as we show next. 
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Constraints on 


Process 


Bound 




r(r — > ISrflV^) 


0.9999 ± 0.0020 




T(r — > v T ev e \ 




T(tc — > nv^) 


1.0017 ± 0.0015 




T(n — > ez/ e ) 


(NM) ee 


r(w -> e p e ) 


0.997 ± 0.010 


(iViVt) Tr 


r(r -> z/ T ez/ e ) 


-1 AAA A 1 A A A £~\ O 

1.0004 ± 0.0023 




-> z^ez/ e ) 


(NN% T 


r(r — > z^ T 7r) 


0.9999 ± 0.0036 




r(vr -> /iz> M ) 


(NN% T 


r(r -> z/ T /iz/ M ) 


1.0002 ±0.0022 


(NN% e 


r(// -> z^ez/ e ) 


(NN% T 

(NM) ee 


-> rz/ T ) 
IW -> ez? e ) 


1.034 ±0.014 



Table 1: Constraints on (NN^) aa from a selection of processes. 
5.4 Rare charged lepton decays 

The leptonic process £ a — > £37 only occurs at one loop, as illustrated in Fig. [31 As the 
photon is on-shell, there are no divergent contributions to the diagram. The 1-loop 
branching ratio in the MUV scheme is given by the same expression [28] as in the 
unitary case, substituting Upmns by N, 



T(4 - ujpup) 32vr (NW) aa (NNi)w ' 

where x k = rnf /Myy with m& being the masses of the light neutrinos and 

10 - 43x + 78a; 2 - 49x 3 ± 4x 4 ± I82 3 In x , . 

S 3(^1)5 • (52) 

Would N be unitary, the x-independent term would vanish exactly through the GIM 
mechanism [2S], for a 7^ (3. With N non-unitary and (NN^) a p 7^ 8 a p, that term 
remains and can be the leading contribution to the branching ratio. With the -very 
accurate- approximation F(x) ~ 10/3, it follows that 



r(4->Vx) 100a \{nn^) 



a/3 1 



12 



(53) 



6 We thank B. Kayser and J. Kersten for reminding us of this departure from the SM prediction 
as a possible signal of non-unitarity and of the importance of radiative corrections in this process, 
respectively. 
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leading to the constraint 



Br ((j, -> ei) < 1.2 ■ 10~ n 
Br{r^ej) < 1.1 • 10 7 , 
Br(r -> /ry) < 4.5 • 10~ 8 , 



(55) 
(56) 
(57) 



\(NN% P \ 2 _ T(4-^7) 967T 



(NN^caiNfifl^p T(4 -> 100a ' 

Strong constraints can now be obtained for the off-diagonal elements of (NN<), using 
the present experimental bounds [30j-[32j 



together with the experimental values Brir — > z/ r /iz/ M ) = 0.1736 ± 0.0006, Br(r — > 
i/ r e!7 e ) = 0.1784 ± 0.0006 and 5r(/i -> v^eV^ « 100% [26]. 

The strong experimental bound on /x — > e7 results i 10 (iViV% (/ , e) < 10" 4 , while 
the other off-diagonal elements are constrained to be less than a few percent. 

Finally, other lepton-flavour violating decays like — > 3£, as well as li — > £j conver- 
sion in nuclei, impose additional constraints close to those above. They may become 
increasingly relevant, depending on the experimental sensitivities attained in the fu- 
ture, as it will be discussed in Sec. [0 



5.5 Summary of constraints on non-unitarity from decays 

All in all, a global fit to the constraints listed in this Section proves that the NN^ 
elements agree with those expected in the unitary case, within a precision better than 
a few percent, at the 90% CL: 

/0.994± 0.005 <7.0-10- 5 < 1.6 - 10 2 \ 
\NN^\ «[ <7.0-10~ 5 0.995 ±0.005 < 1.0 • 10~ 2 . (58) 
\<1.6-10- 2 < 1.0 -10- 2 0.995 ±0.005/ 

7 This strong bound also rules out the possibility of explaining the LSND anomaly with the "zero- 
distance" effect, at least in our minimal scheme. 
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Notice that there is a 2a departure from the value of 1 in the diagonal elements of this 
matrix, stemming from the similar discrepancy in the invisible width of the Z boson, 
not (yet) significant enough to be interpreted SIS db SI gnal of new physics. 

In contrast, there is no direct information from decays on the product N^N, except 
that resulting from the invisible decay width of the Z boson. We can infer strong limits 
on the N^N elements, though, from those in Eq. (1581) . Parametrize the matrix iV as 
N = HV, where V is a unitary matrix and H Hermitian, 

jVjV f = H 2 = 1 + e (59) 

with e = and 

N^N = V ] H 2 V = 1 + V ] eV = 1 + e' . (60) 

e (e') parameterizes the allowed deviation of NN^ (N^N) from the unit matrix. It 
follows that 

i4l 2 s£i4i 2 = £M 2 . (en 

ij a/3 

where the unitarity of V has been used, resulting in the constraint 

i4i<(Ei^i 2 ) 1/2 = - 032 - ( 62 ) 

N^N is thus constrained as follows: 

/1.00± 0.032 < 0.032 < 0.032 \ 

\N*N\ ^ < 0.032 1.00 ±0.032 < 0.032 . (63) 
\ < 0.032 < 0.032 1.00 ±0.032/ 

The results in Eqs. ( 158]) and ( 1631) prove that, within the MUV scheme, unitarity 
in the lepton sector is experimentally confirmed from data on weak decays with a 
precision better than 5%, and does not need to be imposed as an assumption, within 
that accuracy. This means as well that the leptonic unitarity triangles (33J, [M] and 
normalization conditions -corresponding to the elements of iVA^ and N^N- are exper- 
imentally checked with a precision of a few % (or much higher, as for instance for the 
\i - e triangle). 

6 The mixing matrix 

The elements of the mixing matrix obtained from the analysis of neutrino oscillation 
experiments, Eq. (H5l) . can now be combined with the unitarity constraints obtained 
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from weak decays in Eqs. (1581) and (1631) . The resulting mixing matrix in the MUV 
scheme is 

/ 0.75 -0.89 0.45 - 0.65 < 0.20 \ 

\N\ = 0.19 -0.55 0.42- 0.74 0.57- 0.82 . (64) 

\ 0.13 - 0.56 0.36- 0.75 0.54- 0.82 / 

All the elements are now significantly constrained to be rather close to those stemming 
from the usual unitary analysis [35] . 

/ 0.79 -0.88 0.47- 0.61 < 0.20 \ 
\U\=i 0.19 -0.52 0.42- 0.73 0.58 - 0.82 . (65) 
\ 0.20 -0.53 0.44- 0.74 0.56 - 0.81 / 

The constraints resulting for the iV el and N e2 elements are somewhat looser than 
their partners in the unitary analysis. This is due to the large uncertainties allowed for 
the values of the coefficients in Eq. (j4ip . together with the fact that, among all data 
available from solar experiments, we have only included in our analysis the SNO ratio 
of charged to neutral current events. Notice also that the elements of the r-row are 
significantly less bounded than in the unitary analysis, their values being inferred only 
indirectly. 

7 Future experiments 

Matrix elements 

In order to measure independently liV^I 2 and |A^ 2 | 2 without relying on indirect infor- 
mation, a Up disappearance experiment sensitive to Am^ (with the oscillations driven 
by Am| 3 averaged out), as the one proposed in Ref. [M] , would be needed, as suggested 
by Eq. (1391) replacing e by \l. This experiment is quite challenging, requiring an intense 
beam of low energy (~ 500 MeV) and a very long baseline (~ 2000 km). 

Future facilities under discussion, searching for CP-violation in appearance chan- 
nels, include Super-Beams [36], /3-Beams [37] and Neutrino Factories The two 
latter ones may study v e — > and u e —* transitions [39] with great precision, 
while Super-Beams may explore their T-conjugate channels. Unlike disappearance 
experiments, measurements at these facilities will thus be sensitive to the phases of 
the matrix elements, which at present remain completely unknown. Neutrino Factory 
beams would also be energetic enough for the v e — > v T and — > v T oscillation channels 
to be accessible |4T)| |4"1~] . The r-row could thus be tested directly and without relying 
on indirect decay information. 

(iViVt) Me 

As regards unitarity bounds, the constraints on non-unitarity from decays are also 
likely to improve. If no positive signal is found for /i — > ej, the bound on its branching 
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ratio is expected to reach 2 ■ 10~ 14 in the near future [12], which can be translated 
into a unitarity constraint (NN^)^ < 2.9 • 10~ 6 . At a Neutrino Factory, the branching 
ratio for fi — > e'j could be further constrained to < 10~ 15 [13], which would result in 
(iVAT^e < 6.4 • HT 7 , at the 90%CL. 

Important improvements are also expected regarding the bounds for /i to e con- 
version in nuclei. This process is more suppressed than /i — > ej, though, due to the 
extra electromagnetic coupling. In a Neutrino Factory, sensitivities down to 1CT 18 
could be achieved [13] which, translated to (NN^)^ ~ 3.2 ■ 1CT 7 , are only a factor 
two stronger than the bound expected from fi — > cy. Similar ultimate sensitivities are 
being discussed as regards the PRISM/PRIME project [H]. 

(NN% T and (iViVt)^ 

On the other hand, the bounds on rare r decays are not likely to improve much without 
a dedicated facility, since Babar and Belle are now limited by the background and an 
increase in statistics would not improve the relevant measurements [45J. 

In contrast, the possibility of detecting v T at a near detector of a Neutrino Factory 
would allow to improve the bounds on (NN^) eT and (NN^)^. We have considered 
an OPERA-like detector, located at a 100m baseline from a Neutrino Factory bearr|§, 
with a total mass of 4kt and the efficiencies and backgrounds considered in Ref. [41] . 
Assuming a conservative 5% systematic error, the present bounds could be improved 
to (NN^) eT < 2.9 • 10- 3 and (NN*)^ < 2.6 • lO" 3 , at the 90%CL. 



8 Summary and conclusions 

The flavour mixing matrix present in leptonic weak currents may be generically non- 
unitary, as a result of new physics, as for instance that responsible for neutrino masses. 
It is important to know up to which point the values of the matrix elements are allowed 
by data to differ from those obtained in the usual unitary analysis, as putative windows 
of new physics. 

Without attaching ourselves to any particular model, we have studied a minimal 
scheme of unitarity violation -MUV-, considering only three light neutrino species and 
with the usual unitary matrix Upmns replaced by the most general non-unitary one. 

We have first clarified the formalism for studying in this scheme neutrino oscillations 
in vacuum and in matter. It results in "zero-distance" flavour-changing transitions, 
while neutrino propagation in matter exhibits exotic couplings and an active role of the 
neutral current couplings, unlike in the standard unitary analysis. The most relevant 
data from neutrino oscillation experiments have then been used to determine, within 

8 This is only an example of the potential of detecting v T near the Neutrino Factory beam. A 
detailed study of whether the performance of an OPERAdike detector can be extrapolated to the 
neutrino luminosities so close to the source would be required, though. 
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the MUV scheme, the elements of the mixing matrix: about half of its elements remain 
unconstrained at this level. 

Data from weak decays cannot determine the values of the elements of the mixing 
matrix. They provide very stringent tests of unitarity, though, which turns out to be 
satisfied at the level of a few percent. The combination of nowadays' data from weak 
decays and neutrino oscillation experiments results, within the MUV scheme, in a set 
of absolute values for the elements of the mixing matrix, which is very close to that 
from the unitary analysis. It is important to pursue these searches with increasing 
precision, as it may well result in a positive signal instead of constraints, heralding new 
physics. 

Future facilities could establish the first signs of non-unitarity (or new physics in 
general). We have discussed their potential to improve unitarity constraints and de- 
termine all matrix elements without assuming unitarity. While the unitarity bounds 
expected from traditional \i — > experiments will remain very powerful, experiments 
at a Neutrino Factory can set also stringent unitarity constraints, allow a direct deter- 
mination of the elements in the r row of the mixing matrix and furthermore be sensitive 
to the phases of the matrix elements. More generally, the plethora of flavour channels 
and precision measurements, offered by the several types of future facilities under dis- 
cussion, will be very significant in the search for non-unitary effects as windows of new 
physics. 

A final warning is pertinent here. The putative new physics may result in modifica- 
tions to other standard couplings than those considered above, as well as in new exotic 
ones. Barring extreme fine-tunings in some hypothetical model, our approach should 
provide, though, the correct order of magnitude of the minimal unitarity violations 
allowed by data. 
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